Dipole of the luminosity distance: a direct measure of H{z) 
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We show that the dipole of the luminosity distance is a useful observational tool which allows us 
to determine the Hubble parameter as a function of redshift H(z). We determine the number of 
supernovae needed to achieve a given precision for H{z) and to distinguish between different models 
for dark energy. We analyse a sample of nearby supernovae and find a dipole consistent with the 
cosmic microwave background at a significance of more than 2a". 
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One of the biggest cosmological surprises in recent 
years was the discovery that the Universe is presently 
undergoing a phase of accelerated expansion The 
reason for this behaviour is still a complete mystery. 

If the Universe is homogeneous and isotropic on large 
scales, all contributions to the cosmological energy mo- 
mentum tensor are characterised by their energy density 
p[z) and pressure P{z) as functions of cosmic redshift z. 
Accelerated expansion requires that overall p + 3P < 
today. This can be achieved by introducing a so-called 
'dark energy' component with a very negative pressure 
in addition to the usual pressureless matter. One of the 
main challenges of observational cosmology is to char- 
acterise the properties of this dark energy. The homo- 
geneous and isotropic aspects of dark energy are com- 
pletely determined by the equation of state parameter 
Wdc{z) = Pde^z) / pdc{z) which links its pressure and en- 
ergy density. The primary goal of observational dark 
energy studies is to measure the function Wde{z)- 

Current experiments probing the dark energy equation 
of state measure luminosity distances to supernovae or 
the angular diameter distance to the last scattering sur- 
face via the cosmic microwave background (CMB) peak 
positions. These distances are linked to Wdc{z) through a 
double integration, which renders them rather insensitive 
to rapid variations of the equation of state. The required 
modelling can lead to strong biases that are difficult to 
detect and quantify 0. A direct measurement of the 
Hubble parameter H{z) would facilitate the derivation of 
Wdc{z) immensely and allow for a more direct comparison 
with model predictions. As an explicit illustration, let us 
consider a flat universe. The Friedmann equations then 
yield the following relation between Wdc{z) and H{z) 



2/3H'{z) + H^{z) 



m{z) 



{l + zf 



(1) 



Apart from H{z) and its derivative only the parameter 
combination flmHo appears, which will be measured by 
the Planck satellite to an accuracy of about 1%. 

It is possible to obtain H(z) by computing the numer- 
ical derivative of the distance data Isl, but the current 



data leads to a very noisy result. In the future, radial 
baryon oscillation measurements should be able to mea- 
sure H{z) directly. Here we propose an alternative, com- 
pletely independent method based on luminosity distance 
measurements . 

In a previous paper we have considered the luminos- 
ity distance d/, as function of the source redshift z and 
direction n. We have shown that not only the direction- 
averaged luminosity distance. 
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but also its directional dependence, dL(z,n) can be of 
cosmological interest. The directional dependence can 
be expanded in terms of spherical harmonics, leading to 
observable multipoles, Ci{z). In this letter we concen- 
trate on the dipole (corresponding to t ~ 1), 
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Here e is a unit vector denoting the direction of the dipole 
and c?^^^ (z) is its amplitude. 

As we have discussed in (see also 0), for z ^ 0.02 
the dipole is dominated by the peculiar velocity of the 
observer for all redshifts. The lensing contribution to Ci 
is of the order of 10~^ while our peculiar motion induces 
a dipole of Ci ~ 10"^ - lO"*^ for z - 2. At high £'s, i.e. 
small scales the lensing contribution dominates for z ~ 1, 
but this does not interfere with the dipole as it averages 
to zero under the integration Neglecting multipoles 
higher than the dipole we can write the full luminosity 
distance as 
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To derive a formula for d)^ (more details are found 
in 01), we use the luminosity distance to a source emit- 
ting photons at conformal time 77 in an unperturbed 



Friedmann universe, d 
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tion of the observer gives rise to a Doppler effect which 
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is the dominant contribution to the dipole, 

di(77,n)=df (7y)[l-(n.vo)] , (5) 

where Vq is our pecuhar velocity. However, conformal 
time 77 is not an observable quantity, but the source red- 
shift, z = z(r\) + 8z is. Here z(r\) = l/a(ry) — 1 is the 
unperturbed redshift. To first order 
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With d^^\z) = (1 + z){t]o — rf), we have 

^df = {\ + z)-^df +n'\z) and 

5z = — (1 + z)(vo • n) + higher multipoles . (7) 

Here Tiiz) = H(z)/{1 + z) is the co-moving Hubble pa- 
rameter. Inserting this in Eq. Q, we obtain 
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Although Vq is in principle a random variable, we can 
measure it directly from the CMB dipole which is due to 
the same motion. Its magnitude is |vo| — (368 ± 2)km/s 
according to the COBE satellite measurements The 
amplitude of the luminosity distance dipole is then 
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and its direction is e — vo/|vo|. The dipole in the su- 
pernova data gives therefore a direct measure of H{z). 

As a first step we want to test whether there is a dipole 
present at all in the supernova distribution, and if its di- 
rection and magnitude is compatible with expression (jS)) 
(see also 0). Supernova data is conventionally quoted 
in terms of magnitudes rather than the luminosity dis- 
tance. The link between magnitudes and the luminosity 
distance is given by 



X = 51ogio(dL/10pc) 
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where M. is the intrinsic magnitude (however, the abso- 
lute magnitude normalisation is degenerate with log(-ffo) 
and is usually marginalised over). We use the low- 
redshift sample of 44 supernovae assembled by the SNLS 
team Q , together with the supernova locations from Q . 
To the given photometric error we add an error for the 
peculiar velocity of the source of 300 km/s and a constant 
dispersion of Am — 0.12. The latter error ensures a rea- 
sonable goodness-of-fit of both the monopole and dipole 
term. We subtract the monopole of m(2:, n) and find the 
best fit value of Vq for the dipole. In Fig. ^ we show the 
angular uncertainty of Vg. The direction is compatible 




FIG. 1: Direction of the luminosity dipole from a low-redshift 
supernovae sample, in a celestial coordinate system centred 
on the CMB dipole (1 and 2g contours). The two directions 
agree well. 



with the CMB dipole at the Icr level. The magnitude 
of the luminosity dipole gives |vo| = 405 ± 192 km/s, in 
good agreement with the CMB dipole value of 368 km/s. 

Fixing the CMB dipole direction and fitting only the 
ampfitude, we obtain |vo| — 358 km/s with x^^^^-^ — 48.2 
whereas Vq = gives = 52.7. The absence of a dipole 
is therefore disfavored at over 2a. 

Let us estimate the accuracy with which we can de- 
termine H{z) from a measurement of N supernovae in a 
redshift bin [z — dz^z -\- dz]. We assume that the mag- 
nitude of each supernova is known with a precision Am 
(independent of z). We consider vq given by the CMB 
dipole and work with the ansatz Q. The error in the 
magnitude translates into an error on the luminosity dis- 
tance, 

6dL{z,n) — ^^^^^^ dL{z,n)STn{z,n) . (11) 
5 

We add the error into our ansatz, setting 

to(z, n) = to'"' (z) + to'-"'"^ (z)(n • e) + (5m(z, n) (12) 
dL(z,n) = d'^^\z)+d'i\z){n-e)+SdL{z,n). {13) 

We assume that different supernovae are uncorrelated, so 
that the variance of the magnitude is given by 

{Sm{z, n)dm(z, n')) = 47r(Am)^(5^(n - n'). (14) 

The error on the dipole can now be computed using 

Sd^l\z) = ^ I dVL^{n-e)5dL{z,n) . (15) 
47r J 

Its variance is 
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As the monopole is much larger than the dipole we have 
neglected the latter in the previous expression and obtain 
our final formula for the variance of the dipole 
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The absolute error on the dipole is therefore compara- 
ble to the error on the monopole and, not surprisingly, 
the relative error is much larger. 
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We will therefore need a large number of supernovae to 
determine H{z) with reasonable accuracy. 

As the Hubble parameter is inversely proportional to 
the dipole, its relative error is simply 
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This formula is valid for any model of dark energy. 
Once we have measured the luminosity distance, we can 
calculate the monopole and the dipole, deduce the Hub- 
ble parameter, and equation 119(1 gives the error on H{z) 
per supernova at that redshift. We plot the relative error 
on H{z) (which is the same as the relative error on the 
dipole amplitude) in Fig. [21 for a flat universe with a cos- 
mological constant and cold dark matter (ACDM) with 
= 0.7. We use two values for the error on the magni- 
tude. Am = 0.1 and Am = 0.15. This is comparable to 
the accuracy of recent supernova surveys like SNLS 01 . 

In the future it may be possible to control system- 
atic errors much better - indeed the very dipole that 
we want to measure is part of the systematic error bud- 
get of current surveys [loj ! Proper assessment of the 
dipole contribution may therefore also help to measure 
the monopole with higher precision. Our assumed value 
for Am is probably pessimistic as most systematic uncer- 
tainties are expected to affect the overall luminosity at a 
given redshift, i.e. the monopole. The dipole which relies 
on the angular distribution of the luminosity should be 
far more resistant to effects like for example evolution of 
the supernova population with redshift. 

Observing N independent supernovae, the mean error 
on the magnitude is reduced to In Fig. Owe plot the 
number of supernovae needed to measure H[z) with an 
accuracy of 30 %. This number scales quadratically with 
the errors; we need to measure 100 times more super- 
novae to decrease the error by a factor of 10 to 3%. On 
the other hand, if we manage to decrease Am by a factor 
of 10 through an improved understanding of supernova 
explosions and better measurements, then we need 100 
times fewer supernovae. 

One crucial question about dark energy is whether it 
does indeed behave as a cosmological constant or not. 




FIG. 2: We show the relative error in H{z) for one supernova, 
as a function of the redshift, in a flat ACDM universe with 
r^A = 0.7 and for two difi'erent values of the intrinsic error 
Am = 0.1 and Am = 0.15. This represents as well the relative 
error in the dipole d'^l\z). 
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FIG. 3: We show the number of supernovae needed to mea- 
sure H{z) with an accuracy of 30 %, in a flat ACDM universe 
with f^A = 0.7, as a function of the redshift and for two dif- 
ferent values of the intrinsic error Am = 0.1 and Am = 0.15. 



Having measured the dipole at different redshifts, it is 
possible to compare directly the measured value of H{z) 
with the one predicted for a flat ACDM universe. If the 
two do not agree, dark energy must be due to a different 
mechanism, like quintessence or a modification of gen- 
eral relativity. In Fig. ^ we plot the number of super- 
novae needed to distinguish the two cases, by demanding 
that the difference \H{z) — H{^cdm{z)\ be larger than the 
error AH{z). For comparison, the relative difference be- 
tween the Hubble parameter in a flat pure CDM universe 
and in a flat ACDM universe with f^A = 0.7 is 10% at 
z = 0.1, 19% at z = 0.2 and 27% at z = 0.3. 
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Our method tests directly the expansion speed of the 
universe at all the redshifts where we measure the lu- 
minosity distance dipole. Any deviation in H{z) from 
theoretical predictions will be immediately detected. If 
we measure only the usual monopole of the luminosity 
distance then a well-localised deviation may easily be 
smeared out and lost by the additional integration. 
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FIG. 4: We show the number of supernovae needed to differ- 
entiate the measured Hubble parameter H{z) from the the- 
oretical one in a flat ACDM universe Hacdm{z), as a func- 
tion of the relative difference -'^(^~-'^acdm(z)I redshifts 
z — 0.1, z = 0.2 and z = 0.3 and intrinsic error Am = 0.1. 



imal. 

In this letter we have discussed a novel method for mea- 
suring directly the expansion history of the Universe. We 
have shown that the dipole of the supernova distribution 
on the sky is proportional to 1/ H{z). It is therefore pos- 
sible to extract directly H{z) from the dipole. This is 
advantageous compared to the monopole of the luminos- 
ity and angular diameter distance which measure only 
the integral over the Hubble parameter. With a present 
data set of 44 low rcdshift supernovae we have measured 
the dipole and it is in good agreement with the CMB. 

We have also discussed the accuracy with which we can 
measure the Hubble parameter, and found that we need 
a large number of supernovae. However, future planned 
surveys are expected to deliver these. Given that most 
surveys concentrate on high-redshift supernovae while 
the dipole is most useful at moderate redshifts, z ~ 0.5, 
it may be preferable to propose a dedicated low-z super- 
nova survey. 

Finally, the dipole is a quantity independent of the 
monopole. Given a survey with a sufficient number of 
supernovae it is possible to measure both. This improves 
the measurement of the dark energy properties and ad- 
ditionally serves as a cross-check for systematic errors. 

We thank Pierre Astier and Bruno Leibundgut for use- 
ful discussions. This work is supported by the Swiss NSF. 



Looking at the figures it is readily apparent that ac- 
curacy is best at low redshift. This is not necessarily 
a drawback, as dark energy is expected to dominate 
at low redshift and so is best observed in this redshift 
range. However, below z ^ 0.04 non-linear effects prob- 
ably become important which may lead to systematic 
effects in the distribution of supernovae. Future baryon 
oscillation surveys will primarily target higher redshifts, 
z ^ 0.5, where they reach maximum sensitivity. The an- 
gular distribution of low-redshift supernovae is therefore 
a complementary probe. Also the number of supernovae 
needed seems quite realistic. Very large surveys which 
plan to measure of the order of 10** to 10^ supernovae 
are presently discussed JTT) . 

As a final remark, even though uniform sky coverage 
is not essential, a survey designed to measure the dipole 
should optimally cover a large part of the sky. If we only 
observe supernovae in one small patch, it may be difficult 
to extract the dipole without contamination from lensing 
which dominates over the dipole for £ ~ 100 and z ^ 1 
(see 0). If possible, the observed supernovae should 
cover the regions of the sky aligned and anti-aligned with 
the CMB dipole where the luminosity difference is max- 
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